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摘 EB: 作物 模型 的 研究 涉及 作物 生长 发 育 的 复杂 过 程 ， 空 间 上 从 分 子 到 细胞 、 组 织 、 器 官 、 个 体 、 和 群体 
等 不 同 尺度 ， 时 间 尺 度 上 可 以 从 秒 到 年 。 基 于 不 同 的 研究 需求 ， 切 换 作物 模型 尺度 ， 可 使 得 作物 模型 的 适 
用 性 更 广泛 灵活 。 其 中 ， 如 何 从 群体 尺度 的 作物 模型 转 和 人 个体 尺度 的 作物 模型 是 本 研究 的 内 容 。 本 研究 基 
于 四 个 玉米 品种 的 两 个 处 理 (灌溉 和 而 养 ) 的 已 有 的 实验 数据 和 基于 这 些 数 据 的 DSSAT 系 统 的 模拟 数据 ， 
校准 功能 结构 模型 GreenLab 的 参数 ， 以 计算 结果 一 致 为 指标 ， 探 索 不 同 空间 尺度 模型 建立 接口 的 方法 ， 比 
较 不 同 模型 的 特点 。 结 果 表 明 ，GreenLab 模型 可 以 复 现 DSSAT 系 统 的 模拟 数据 和 实际 测量 数据 ， 进 一 步 可 
以 反 演出 各 种 器 官 之 间 生 物 量 的 分 配 并 进行 三 维 可 视 化 展示 。 最 后 讨论 了 不 同 空间 尺度 模型 结合 的 优势 及 


应 用 领域 。 
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English abstract) 
1 引言 展 以 及 植物 生长 机 理 的 深入 研究 ， 学 科 间 的 交叉 


促进 了 作物 模型 的 发 展 。 利 用 计算 机 对 作物 生长 

模型 是 测试 假设 、 综 合 知识 、 描 述 和 理解 复 发育 过 程 及 其 与 环境 的 动态 关系 进行 系统 的 分 析 

杂 系 统 并 比较 不 同方 案 的 有 力 工具 ， 可 用 于 决策 。 和 定量 的 描述 ， 可 预测 作物 的 阶段 发 育 、 形 态 发 

支持 系统 、 温 室 气 候 控 制 和 预测 以 及 生产 计 干 物质 积 分 配 及 产量 ， 通过 综合 知识 的 表 
划 ”…。20 世 纪 60 年 代 以 来 ， 随 着 信息 技术 的 发 ， 深 入 洞察 作物 生长 过 程 OU 
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作物 模型 的 研究 涉及 作物 不 同时 空 尺 度 的 复 
杂 生 长 过 程 ， 空 间 尺 度 从 分 子 到 细胞 、 组 织 、 带 
官 、 个 体 、 群 体 等 ， 时 间 尺 度 包 括 秒 、 天 、 年 
等 中。 作物 模型 尺度 的 选择 取决 于 研究 的 对 象 、 
过 程 和 人 类 认 知 的 程度 等 。 如 图 1 所 示 为 作物 模 
型 从 基因 到 生态 系统 的 时 空 尺度 的 划分 ”“， 有 助 
于 清晰 定位 所 研究 的 模型 尺度 ， 保 证 在 同一 尺度 
下 讨论 问题 。 不 同 尺 度 模型 的 融合 可 以 取 长 补 
短 ， 使 得 模型 的 适用 性 更 广泛 。 


10° 10° 10° 10° 108 10" 10° 10" 
时 间 /s 


图 1 作物 模型 研究 在 不 同时 间 和 空间 尺度 的 划分 
Fig. 1 Scales of time and space for different crop models 

作物 模型 可 分 为 描述 性 模型 (Descriptive 
Model) 和 人 解释 性 模型 (Explanatory Model)。 前 
者 又 称 为 统计 模型 、 回 归 模 型 、 经 验 模型 或 黑箱 
模型 ， 直 接 建立 输入 与 输出 之 间 的 关系 ， 不 关心 
内 在 形成 机 理 。 后 者 又 称 为 基于 过 程 的 模型 
(Process-Based Model) ， 是 指 能 定量 地 和 动态 地 
描述 作物 生长 、 发 育 和 产量 形成 过 程 及 其 对 环境 
反应 的 计算 机 模拟 程序 “， 通 常 是 在 群体 水 平 对 
作物 生长 的 描述 。 和 常 见 的 基于 过 程 的 作物 模型 ， 
如 CERES " , TomSim " , STICS (Simulateur 
mulTldisciplinaire pour les Cultures Standard) ”和 
APSIM (Agricultural Production Systems sIMula- 
tor) " 等， 模拟 作物 的 光合 作用 过 程 、 同 化 物 的 
分 配 过 程 等 ， 主 要 通过 叶 面 积 指数 (Leaf Area 


难以 描述 作物 个 体内 结构 的 变化 及 其 对 产量 的 影 
啊 。 以 作物 模型 为 核心 的 农业 技术 转移 文 持 决策 
系统 (Decision Support System for Agrotechnolo- 
gy Transfer, DSSAT) 以 天 、 日 为 步 长 模拟 作物 
生长 ， 按 照 作物 生育 期 详细 描述 作物 生长 发 育 过 
fe, 包括 发 芽 到 开花 、 叶 片 出 现 、 开 花 时 期 、 籽 
粒 灌浆 、 生 理 成 熟 和 收获 等 ， 并 且 可 响应 许多 因 
素 ， 包 括 作物 遗传 特性 、 管 理 措施 、 环 境 、 氮 素 
和 水 分 的 胁迫 、 病 虫害 等 ， 主 要 用 于 农业 试验 分 
析 、 农 业 产量 预报 、 农 业 生产 风险 评估 、 和 气候 对 
农业 的 影响 评价 等 。 

由 于 作物 结构 的 重要 性 ,衍生 了 个 体 尺 度 的 
植物 功能 结构 模型 (Functional-Structural Plant 
Model, FSPM) 来 模拟 作物 形态 结构 、 生 物 量 
的 产生 和 分 配 以 及 两 者 的 内 在 联系 。FSPM 是 在 
单个 器 官 到 植物 个 体 层次 对 植物 这 一 复杂 系统 的 
建 模 与 仿真 ""， 在 过 去 的 二 十 多 年 间 得 到 了 迅 
速 的 发 展 ， 成 为 植物 建 模 的 研究 热点 之 一 。 

GreenLab 模 型 是 作物 功能 结构 模型 之 一 ， 考 
虑 同类 器 官 顺序 产生 的 特点 ， 并 采用 公共 池 和 源 
库 关系 的 概念 ， 既 在 叶 元 尺度 上 (器官 ) 对 植物 
的 生长 进行 建 模 ， 同 时 又 保持 了 与 基于 过 程 的 作 
物 模型 (群体) 的 兼容 性 。GreenLab 模型 基于 植 
物 自 动机 (Automaton) 2 ， 通 过 递归 算法 计算 
每 个 时 间 步 长 产生 的 器 官 数量 ， 再 根据 相对 库 强 
(Relative Sink Strength) 和 器 官 数量 ， 将 从 公共 
池 中 获取 的 生物 量 分 配 到 不 同类 别 的 需 官 中 。 该 
过 程 通过 数学 公式 来 撕 述 ， 不 需要 逐个 器 官 地 来 
模拟 生物 量 的 分 配 ， 因 此 计算 速度 快 ， 所 需 时 间 
少 ， 这 是 GreenLab 模型 的 优势 之 一 。 此 外 ， 
GreenLab 模 型 的 特色 在 于 通过 实测 植物 各 个 器 官 
的 重量 等 数据 可 反 求 影响 生物 量 产生 和 分 配 的 模 
型 源 库 参 数 , CHM RK pe W 
AUS. BOI HO 7 等 十 几 种 作物 进行 模型 


Index, LAI) 来 预测 每 平方 米 的 生物 量 产 量 ， 同 
时 考虑 辐射 、 温 度 和 基本 投入 《〈 灌 溉 和 化 肥 等 ) 
的 影响 。 这 类 模型 大 多 将 不 同类 型 的 器 官 作 为 一 
个 整体 考虑 ， 关 注 作 物 的 总 果实 重量 、 叶 重 等 ， 


应 用 人 研究。 但 是 ，GreenLab 模 型 将 环境 因素 
的 影响 简化 为 一 个 环境 因子 E 来 描述 ， 不 能 很 好 
地 模拟 气息、 土壤 以 及 管理 措施 等 对 作物 产量 的 
作用 。 


ChinaXiv 合 作 期 刊 


79 


202302.00207v1 


chinaXiv 


近年 来 ， 随 着 作物 表 型 信息 获取 设备 的 迅速 
发 展 ， 已 实现 采集 作物 株 高 、 叶 面积 、 籽 粒 数 ， 
以 及 生理 和 光合 等 各 种 表 型 信息 。 因 此 ， 如 何 对 
获取 的 大 量 表 型 数据 进行 分 析 是 当前 正 待 解决 的 
问题 "”。 不 同 作物 模型 的 融合 ， 有 助 于 更 深入 
地 理解 作物 行为 ， 更 大 程度 发 挥 模型 作用 。 本 人 研 
究 采 用 DSSAT 系统 与 GreenLab 模型 进行 群体 尺 
度 和 个 体 尺 度 模型 的 融合 研究 ， 探 索 以 不 同 空间 
尺度 和 方法 所 构建 的 模型 进行 交互 的 可 能 性 ， 使 
得 既 可 以 模拟 作物 遗传 特性 、 管 理 措施 、 环 境 等 
的 影响 ， 同 时 也 可 以 模拟 植物 个 体内 各 器 官 生长 
和 发 育 的 动态 过 程 ， 根 据 对 植物 的 更 深层 次 的 了 
解 来 详细 描述 产量 的 组 成 ， 并 对 作物 进行 三 维 展 
示 ， 与 表 型 信息 做 对 比 。 通 过 两 类 模型 的 融合 ， 
为 植物 育种 、 表 型 鉴定 和 作物 系统 优化 提供 
EF. 


2 材料 与 方法 


本 研究 总 体 思 路 如 图 2 所 示 。 模 型 融合 的 思 
想 考 虑 到 发 展 模 型 的 人 员 、 模 型 实现 方法 的 不 
同 ， 且 自 成 体系 ， 选 择 以 统一 的 输出 结果 进行 连 
接 ， 而 不 是 尝试 在 模型 构成 模块 进行 融合 。 
先 ， 利 用 DSSAT 系统 输出 模拟 的 每 日 玉米 的 叶 
重 、 茎 重 和 玉米 棒 重 等 ;再 基于 这 些 数据 ， 结 合 
GreenLab 模型 已 有 的 玉米 的 结构 信息 与 库 强 变异 
参数 等 结果 ， 计 算 模 型 的 源 库 参数 ， 输 出 不 同 叶 
位 、 不 同 生 长 阶段 的 各 器 官 的 详细 数据 ， 以 及 玉 
米 植株 的 三 维 结构 图 。 


2.1 数据 选取 与 处 理 


2.1.1 试验 数据 选取 

试验 数据 包括 两 部 分 。(1) Soler ^ ^" 2002 
年 在 巴西 皮 拉 西 卡巴 ( 西 经 22.7?， 南 纬 47.4°， 
iik 580mm) 进行 的 四 个 玉米 品种 (AG9010、 
DKB333B、DASCO32 和 EXCELER) 两 个 处 理 ， 
BWA (Irrigated, I) 与 雨 养 (Rained, R) 的 
田间 试验 测量 数据 ， 包 括 土壤 、 气 候 、 管 理 措施 
等 ， 以 及 七 次 定期 (3 月 12 日 、3 月 21 日 、3 月 


历史 天 气 、 土 壤 、 作 物品 种 、 管 理 措施 
输入 
Y 
DSSAT 系 统 
输出 
v 
FAE, SB, GEOK Ex 
输入 
x 
GreenLab 模 型 


-— < 


模型 参数 计算 


Y Y 
不 同位 置 、 不 同 生长 阶段 的 每 个 节 x j às 
MIHE ZE, ORR 玉米 的 三 维 结构 展示 


图 2 DSSAT 4 GreenLab 融合 技术 路 线 图 
Fig. 2 Framework ofthe integration between DSSAT and 


GreenLab 


27 日 、4 月 7 日 、4 月 18 日 、 4 月 28 日 和 5 月 14 
日 ) 测量 的 玉米 LAI 和 地 上 部 总 重量 ， 共 包含 8 
个 对 比 实验 。 其 中 ，AG9010 为 短 季 品种 (从 种 
植 到 吐 丝 期 为 904 生 长 度 日 (Growing Degree- 
Day，GDD) ， 基 温 为 8"C) ，DASCO32 与 EX- 
CELER 为 两 个 短 季 品种 (995 GDD), DKB333B 
为 常规 品种 (1037 GDD) 2 2 。 所 有 试验 为 4 个 
重复 的 随机 区 组 设计 。 每 个 地 块 长 20 m, 行距 
0.8 m， 共 四 行 ， 种 植 密度 为 5 k/m; (2) 基于 
(1) 中 的 数据 ， 利 用 DSSAT 系 统 模拟 4 个 玉米 品 
种 的 每 日 数据 ， 包 括 每 日 叶 、 节 、 玉 米 棒 和 地 上 
部 总 重 等 。 本 研究 选取 了 叶片 数 分 别 为 7、10、 
15、20、23、27、30、33、36、39 和 40 共 11 个 
生长 阶段 的 DSSAT 系统 模拟 数据 用 于 GreenLab 
模型 的 模拟 计算 。 

基于 (2) 中 的 DSSAT 系统 模拟 的 每 日 数据 
拟 合 GreenLab 模型 的 参数 ， 模 拟 LAI、 玉 米 叶 、 
节 和 玉米 棒 的 重量 ， 以 及 地 上 部 总 重量 ， 并 比较 
四 个 玉米 品种 和 两 个 处 理 的 DSSAT 系统 和 
GreenLab 模 型 模拟 的 数据 与 实验 测量 数据 之 间 的 
差异 性 ， 分 析 模 型 结合 的 可 行 性 。 
2.1.2 ”数据 处 理 

DSSAT 模 型 以 日 为 单位 进行 模拟 ， 而 Green- 
Lab 模 型 以 叶 元 为 单位 进行 模拟 ， 每 长 出 一 片 叶 
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子 即 被 认为 是 一 个 生长 周期 (Growth Cycle, 
GC)， 即 出 叶 间 隔 (phyllochron) ， 进 行 模 拟 。 
因此 ， 需 要 将 日 转换 为 生长 周期 。 通 过 将 日 历时 
间 CX) 转换 为 积温 (温度 总 和 ) 来 调节 植物 的 
发 育 ， 从 而 可 以 定义 周期 的 概念 。 

在 GreenLab 模型 中 ， 一 个 生长 周期 中 一 个 
分 生 组 织 产 生 一 个 叶 元 ， 叶 子 接收 与 该 周期 中 累 
积 的 光 辐 射 成 比例 的 能 量 。 根 据 DSSAT 中 的 每 
日 数据 与 积温 ， 计 算出 GreenLab 模型 中 对 应 的 
生长 周期 。 根 据 种 植 密度 计算 出 每 株 植 物 各 类 器 
官 的 重量 。 

为 同步 植物 功能 结构 模型 中 的 生长 和 发 育 ， 
GreenLab 中 的 生长 模型 采用 与 发 育 模型 相同 的 时 
间 尺 度 ， 即 生长 周期 ， 从 而 简化 了 模型 的 计算 。 
在 每 个 周期 中 ,分 生 组 织 产 生 新 的 叶 元 CES); 
同时 ， 叶 片 进行 光合 作用 产生 生物 量 并 根据 各 自 
的 库 强 分 配 到 植株 的 各 个 带 官 中 (生长 )。 因 此 ， 
基于 同一 时 间 尺 度 模拟 植物 的 发 育 和 生长 。 这 种 
方法 使 得 当 植 物 的 顶端 发 育 停 止 〈 雄 穗 出 现 ) 
时 ， 仍 然 可 以 确定 其 时 间 斥 度 。 生 长 周期 数 与 积 
温 之 间 的 线性 关系 是 在 植物 发 育 阶 段 建立 的 ， 然 
后 外 推 到 植物 生长 阶段 。 

此 外 ，DSSAT 模 型 中 的 重量 单位 为 gm， 需 
要 根据 种 植 密度 计算 出 每 株 植物 的 各 类 需 官 的 重 
量 ， 本 研究 主要 考虑 叶 、 节 和 果 的 重量 。 


2.2 模型 参数 设置 


基于 过 程 的 作物 模型 原理 是 在 日 历 (天 ) 的 
每 个 时 间 步 长 计算 单位 面积 土地 截获 的 光 所 产生 
的 生物 量 。 植 物 吸收 的 光 辐 射 量 可 以 根据 Beer- 
Lambert 定 律 计算 ， 公式 为 : 

II, = ee (1) 

其 中 ，k 为 消光 系数 ; 7 为 冠 层 深度 为 工时 的 
辐射 水 平 ，W/m ; 下 分别 为 冠 层 最 上 方 的 辐射 水 
Æ, Wim’; k 值 一 般 在 0.5 至 0.8 的 范围 内 路。 光 
吸收 量 随 叶片 面积 的 增加 而 增加 ,但 叶片 遮挡 会 
减少 光 的 截获 。 

生物 量 可 通过 公式 (2) 计算 。 


Q(t) = RUE*(1 - exp(-k-LAI (t))) (2) 

其 中 ，O (t) 为 在 时 间 t 产 生 的 生物 量 ，g; 
RUE 表示 光 能 利用 效率 (Radiation Use Efficien- 
cy); LAI (t) 是 指 在 时 间 t 测 量 的 群体 叶 面 积 指 
数 。 通 过 对 时 间 进 行 求 和 ， 可 以 通过 累加 计算 总 
«Ey. 

GreenLab 模型 是 基于 生物 量 公 共 池 7 以 及 
源 库 关系 假设 的 模型 ， 并 使 用 了 作物 生长 年 龄 和 
生理 年 龄 的 概念 。 作 物 生长 年 龄 指 植 物 所 经 历 的 
生长 周期 数 ， 所 有 器 官 都 有 生长 年 龄 。 作 物 生理 
年 龄 用 来 区 分 植物 中 不 同类 型 的 枝 ， 主 葵 生 理 年 
龄 为 1， 分 枝 生 理 年 龄 为 2。 

GreenLab 模 型 中 的 时 间 步 长 转化 为 生长 周期 
表示 。 计 算 生 物 量 生产 的 公式 为 : 


E'S, k max (t, £,) me argc; 
Q(t)= r top AS sere (i): 
+ Pi(j-i+1)Q(j-1) 
3 
2 DG) | à 


其 中 ,EE 为 潜在 光合 生产 能 力 ， 即 将 环境 因 
素 的 影响 简化 为 一 个 参数 来 表示 ; 1 为 植物 的 生 
长 年 龄 ， 即 植株 的 生长 周期 ，GC; o HERE 
龄 ， 从 1 到 mxg， 玉 米 为 单 蕉 植物; i, 7 为 不 同 
器 官 的 生长 年 龄 ，GC; zt 为 叶子 的 功能 周期 ， 可 
根据 经 验 或 观测 值 设 定 ; 5, 为 植株 投影 面积 ， 
m”， 取 值 为 种 植 密度 的 倒数 "站 ; s 为 叶片 厚度 ， 
一 般 为 常数 ， 可 通过 测量 的 叶片 重量 除 以 叶 面 积 
计算 gem’; 7 为 光 利 用 系数 ; ky Beer-Lambert 
定律 的 阻力 系数 ,一般 为 0.7'，”; NOO 为 周期 i 
时 的 叶子 数 ， 用 于 计算 植株 总 的 叶片 面积 ; 
P,(j -i+ 1) 为 功能 叶片 的 库 强 ; D G) 为 周期 / 
的 植物 总 需求 ; 0(j - 1) 为 上 一 周期 产生 的 植物 
总 生物 量 ，g; O (0) 为 种 子 生 物 量 ，g。 

为 了 计算 植物 的 叶 面积 ， 需 要 将 功能 性 叶子 
的 重量 相 加 ， 再 除 以 它们 的 厚度 e《〈 假 定 为 常 
数 )。 可 以 通过 公式 (4) 计算 功能 叶 面 积 
Sf (Ws 


1 mx max (t, t.) : B 
Sf()- c7 Eua NO aly) (4) 
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其 中 ， /为 植物 生长 年 龄 ，GC; i 为 叶片 生 
长 年 龄 ，GC; N? (i, 1) 为 植物 生长 年 龄 为 :， 叶 
片 生长 年 龄 为 i 时 叶片 的 数量 , 个; gr (i, t) AH 
物 生 长 年 龄 为 叶片 生长 年 龄 为 ;时 叶片 的 生物 
量 ，g。 在 生长 年 龄 为 ;的 一 组 叶片 中 ， 记 录 了 叶 
片 的 数量 Ns (i,t) 及 其 各 自 的 生物 量 qs(i,1)。 根 
据 植 物 的 年 龄 是否 超 过 叶片 1 的 寿命 ， 需 要 划 
分 为 两 种 情况 来 计算 面积 。 实 际 上 ， 只 要 没有 叶 
子 死亡 ，!< 大 的 情况 就 相当 于 植物 生长 的 开始 。 
在 男 一 种 情况 1>t 下 ， 只 有 生长 年 龄 小 于 1 一 t+ 
1 的 叶子 仍然 起 作用 ， 这 部 分 通过 公式 (4) 中 的 
第 二 个 和 的 索引 计算 。 

由 于 玉米 生长 后 期 叶 和 节 的 生物 量 会 转移 到 
果实 中 ， 因 此 ， 还 需要 考虑 物质 再 分 配 过 程 
(Remobilization) ， 即 成 熟 的 需 官 (HUN) E 
新 分 配 部 分 生物 量 给 果 ， 计 算 公 式 如 下 。 

g(t) = au [17 kmat (1 = Kio) 7) (5) 

其 中 ，1 为 植物 的 生长 年 龄 ，GC; qx 
顺 官 的 最 大 重量 ，g; ,为 该 带 官 转移 的 生物 量 
的 比例 ，%， 可 通过 测量 数据 进行 估计 ; ,为 生 
物 量 转 移 速 度 ，g/GC， 由 模型 反 求 得 到 ; o 表 示 
器 官 类 型 ， 本 研究 中 为 叶 b AI io 
2.2.1 模型 生长 参数 设置 

根据 种 植 密度 4 为 5 株 /mr 设 定 5 的 初始 值 
为 Wd ~ 2000， 实 际 根据 拟 合 结果 进行 调整 。 根 
据 DSSAT 中 的 光 能 利用 效率 (RUE: 4.2 g/MJ) 
以 及 GreenLab 模 型 中 的 经 验 值 设 定 与 DSSAT 中 
相对 应 的 r 值 为 50; 根据 DSSAT 模 拟 的 数据 ， 设 
定 叶 、 节 和 果 的 扩展 周期 分 别 为 15、25 和 20, 
功能 周期 都 为 25， 以 及 用 于 计算 库 强 变异 函数 的 
ZU. 设 定 叶 片 的 相对 库 强 为 1。 

222 ”发 育 参 数 设 置 

对 于 玉米 ， 每 个 叶 元 包含 一 个 叶子 、 一 个 节 
或 果 。 根 据 DSSAT 中 的 数据 ， 最 终 的 叶子 数 为 
19~-21 个 ; 根据 DSSAT 模型 中 的 数据 ， 四 个 玉 
米 品种 都 为 一 个 果 的 类 型 ， 位 于 主 葵 的 第 16 或 
17 个 节 上 。 


2.3 模型 参数 拟 合 


根据 公式 〈3)， 将 其 中 的 > 值 固 定 ， 拟 合 环 
境 参数 EE 和 库 强 参数 P，(o 为 叶子 b、 节 i 和 果 了 )， 
5 的 值 根据 模拟 结果 进行 适当 调整 。 根 据 经 验 ， 
环境 的 变化 对 GreenLab 模 型 中 的 库 强 参数 影响 
不 大 '""。 因 此 ， 尝 试 采用 一 个 环境 参数 的 平均 
值 E£， 而 不 是 读 取 环 境 序 列 进行 模拟 。 此 外 ， 还 
需要 拟 合 物质 再 分 配 (Remobilization) 的 两 个 
BRK, AM ko FOP, PEE IK, Ak, HRA SE Soe 
数据 估算 得 到 ， 均 为 0.8; 转移 速度 参数 ,和 
由 模型 拟 合 得 到 。 


3 试验 结果 与 分 析 


3.1 试验 结果 


3.1.1. 株 高 和 单位 面积 最 终 产量 

四 个 玉米 品种 的 两 个 处 理 收获 时 的 株 高 和 单 
位 面积 最 终 产量 数据 如 表 1 所 示 。 灌 小 处理 的 株 
高 和 单位 面积 最 终 产量 均 高 于 雨 养 处 理 的 植株 。 
其 中 ， 品 种 EXCELER 两 个 处 理 的 产量 均 最 大 ， 
品种 DASCO32 次 之 ， 品 种 AG9010 最 小 。 绪 果 
表明 ， 灌 溉 可 以 显著 提高 四 个 玉米 品种 的 产量 ， 
且 EXCELER 是 最 优 品种 。 
表 1 四 个 玉米 品种 两 个 处 理 的 株 高 和 单位 面积 最 终 产 量 
Table 1 Plant heights and final yield of the two treatments 


for the four maize cultivars 


i 株 高 /cm 产量 /gm-) 
icai AK me GE UPS 
AG9010 200 175 498.6 404.4 
DASCO32 217 206 504.7 410.9 
EXCELER 230 227 530.6 485.9 
DKB333B 213 202 513.9 382.2 


3.1.2 DSSAT 模 拟 的 各 类 器 官 重量 

四 个 玉米 品种 的 两 个 处 理 收获 时 的 叶子 、 
节 、 果 各 类 器 官 重量 如 图 3 所 示 。 洪 溉 处 理 的 各 
类 器 官 重量 均 高 于 不 灌溉 处 理 。 但 DSSAT 模拟 
的 DKB333B 的 产量 最 大 ， 与 试验 测量 结果 偏差 
较 大 。DASCO32 和 EXCELER 的 产量 相近 ， 
AG9010 的 产量 最 小 ， 与 表 1 试 验 结果 一 致 。 
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x v —Q" unu. 一 一 AG9010 计算 值 € ”AG9010 测量 值 
EZZA "x MN 4H JOKHGE NENN 6 DOE 一 一 DASCO32 计算 值 e DASCO32 EET 
300 一 一 EXCELER 计算 值 “@ — EXCELER 测量 值 
一 一 DKB333B 计算 值 €  DKB333B 测量 值 
50 
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40 + 
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M 
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AGR EXCR DASR DKBR AGI EXCI DASI DKBI 
品种 和 处 理 
i: AG, EXC, DAS, DKB 为 玉米 的 四 个 品种 ; R 为 雨 养 , LAB 
图 3 DSSAT 系统 模拟 的 四 个 品种 两 个 处 理 的 玉米 收获 阶 
段 各 类 器 官 ( 叶 、 节 、 果 ) 的 重量 及 总 重量 
Fig. 3 Weights of organs (leaves, internode, fruits) and total 


weight for the four maize cultivars and two treatments 
3.2 模型 模拟 结果 


3.2.1 两 个 模型 时 间 尺 度 的 转换 

四 个 玉米 品种 灌溉 和 两 养 处 理 的 叶子 数 和 理 
论 生长 周期 数 随 积温 的 变化 相同 (图 4)。 因 此 ， 
只 给 出 了 灌溉 处 理 的 四 个 玉米 品种 的 生长 周期 随 
积温 的 变化 。 不 同 品种 间 的 叶子 数 无 显著 差异 
(19~21)， 且 三 个 短 季 品种 AG9010、DAS CO32 
45 EXCELER 所 需要 的 生长 周期 数 无 差异 ， 仅 党 
规 品 种 DKB 333B 的 生长 周期 数 略 长 ， 但 差异 不 
显著 (F=0.097, P>0.5, df=3). 

根据 温度 和 叶片 数量 可 以 计算 出 叶 间 隔 。 对 
于 玉米 ， 积 温 与 所 发 育 的 叶 元 数量 之 间 呈 线性 关 
系 。 因 此 ， 可 以 计算 出 以 积温 表示 的 生长 周期 ， 
以 及 各 生长 周期 对 应 的 叶 元 数量 。 

当 叶 片 数量 不 再 增加 后 ， 玉 米 到 达成 熟 期 ， 
该 阶段 虽然 不 再 产生 新 的 叶 元 ， 但 生长 仍 在 继 
续 。 可 以 根据 从 发 芽 到 植物 项 端 抽穗 阶段 的 积温 
和 叶片 数量 之 间 的 二 元 线性 关系 (斜率) 来 确定 
理论 的 生长 周期 数 。 即 将 时 间 离 散 化 ， 根 据 发 育 
阶段 确定 的 线性 关系 来 计算 生长 阶段 的 理论 生长 
周期 数 (图 4)。 
3.2.2 ”模型 源 库 参 数 

对 于 每 个 实验 ， 通 过 将 参数 %、r， 以 及 库 


0 500 1000 1500 2000 
积温 /("C-d) 


图 4 四 个 玉米 品种 灌溉 处 理 的 生长 周期 和 叶片 数 随 积温 
的 变化 
Fig. 4 Growth cycle and number of leaves with thermal time 


for the four maize cultivars and irrigation treatments 


强 变 异 参数 设置 为 固定 值 ， 可 以 拟 合 环境 参数 
巨 、 库 强 参 数 (PAM PO 以 及 物质 再 分 配 参 数 ， 
All ky 

拟 合 的 参数 值 如 表 2 所 示 ， 可 以 看 出 四 个 玉 
米 品 种 灌溉 处 理 的 已 、 已 、 已 和 8, 的 值 均 大 于 不 
灌溉 处 理 。 其 中 ，EXCELER 两 个 处 理 的 果实 
库 强 (Pj) 最 大 ，DASCO32 和 DKB333B 次 
之 ，AG9010 最 小 。 该 结果 与 试验 结果 一 致 ( 见 
表 1)， 四 个 玉米 品种 两 个 处 理 的 模拟 值 与 测量 值 
的 相关 系数 均 为 0.99， 表 明 模 型 的 模拟 精度 


表 2 四 个 玉米 品种 两 个 处 理 的 拟 合 参数 值 
Table 2 Parameter values of the two treatments for the 


four maize cultivars 


品种 ”处理 E P Pp ky A/(gGC?) Sya 


灌溉 0.17 0.96 243.3 0.014 0.030 3000 


AG9010  ., 

NH 0.15 0.76 71.9 0.016 0.035 2500 

灌溉 0.23 1.26 692.6 0.014 . 0.030 2200 
DASCO32 __ 

HA 0.17 0.89 112.9 0.013 0.033 2150 

灌溉 0.201 1.29 865.8 0.015 0.035 2500 
EXCELER . |. 

A 0.17 1.03 143.2 0.015 0.050 2200 

灌溉 0.21 1.22 577.6 0.05 0.030 2100 
DKB333B __.. 

H 0.18 0.94 114.9 0.016 0.032 2000 


3.3 两 个 模型 模拟 结果 与 测量 结果 比较 


由 于 田间 实验 只 测量 了 LAI 和 地 上 部 总 重 
量 ， 因 此 ， 将 两 个 模型 的 LAI 和 地 上 部 总 重量 的 
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模拟 结果 与 实验 数据 进行 对 比 ， 以 验证 模型 的 准 
确 性 。 
3.3.1 LAI 的 比较 

图 5 为 不 同 阶段 LAI 的 测量 值 与 DSSAT 和 
GreenLab 模拟 值 的 比较 。 由 图 可 见 GreenLab 
模型 计算 的 LAI 与 实际 测量 值 的 一 致 性 更 好 ， 且 
AG9010 和 DKB 333B 两 个 品种 的 DSSAT 模拟 值 
在 后 期 偏 高 。 结 果 表 明 ， 基 于 DSSAT 的 模拟 数 
据 ， 通 过 GreenLab 模型 可 以 更 好 地 模拟 
LAI 值 。 


R Cmp — R Sim e RObs e I Obs ——1Cmp —— I Sim 


AG9010 


DKB333B 


4 
DASCO32 EXCELER 


0 10 20 30 40 0 10 20 30 40 
生长 周期 /GC 


i£; Cmp, Sim fe Obs 分 别 表 示 GreenLab 模 拟 值 ,DSSAT 模 拟 值 
和 测量 值 iI 和 及 分 别 表示 灌流 和 雨 养 处 理 
图 5 四 个 玉米 品种 两 个 处 理 的 LAI 值 
Fig.5 LAI ofcomputed for the four maize cultivars and 


two treatments 


3.3.2 ”地 上 部 总 重量 的 比较 

图 6 为 不 同 阶段 的 地 上 部 总 重量 的 测量 值 与 
DSSAT 和 GreenLab 模 拟 值 比较 。 二 者 的 灌溉 和 
雨 养 两 个 处 理 的 模拟 结果 都 与 实验 测量 的 地 上 部 
总 重量 均 具 有 很 好 的 一 致 性 。 结 果 表 明 ， 两 种 模 
型 可 以 很 好 地 模拟 这 些 值 ， 具 有 一 致 的 输出 。 


3.4 GreenLab 与 DSSAT 模 拟 结果 比较 


3.4.1 不 同 生 长 阶段 各 类 器 官 重量 模拟 
GreenLab 模 拟 的 四 个 玉米 品种 两 个 处 理 的 不 
同 生 长 阶段 各 类 需 官 重量 结果 均 非 常 好 ， 所 有 模 


——R.Cmp 一 -RSsim e Robs e I Obs —— I Cmp — I Sim 
300 300 
AG9010 DKB333B 


Top Weight/g 


EXECLER 


DASCO32 


0 10 20 30 40 0 10 20 30 40 
生长 周期 /GC 


注 : Cmp,Sim 和 Obs 分 别 表 示 GreenLab 模拟 值 ,DSSAT 模拟 值 和 
测量 值 ;I 和 及 分 别 表 示 灌 溉 和 雨 养 处 理 
图 6 四 个 玉米 品种 两 个 处 理 的 地 上 部 总 重量 


Fig. 6 Top weights of computed for the four maize cultivars 


and two treatments 


拟 的 尼 均 为 0.99， 如 岁 7 所 示 。 结 果 表 明 Green- 
Lab 模型 可 以 很 好 地 计算 这 些 值 ， 计 算 结果 与 
DSSAT 的 模拟 结果 一 致 。 
3.4.2 各 类 器 官 重量 随 位 置 和 生长 阶段 的 变化 
模拟 

GreenLab 模型 是 基于 需 官 尺度 的 模型 ， 可 以 
模拟 各 类 器 官 重 量 随 位 置 (图 8 (a) ) 和 生长 阶 
Be (图 8 (b) ) 的 变化 情况 。 叶 子 和 节 的 重量 随 
叶 位 和 生长 阶段 均 呈 先 增加 后 减 小 的 趋势 。 
3.4.8 ”生物 量 生产 和 需求 的 变化 模拟 

GreenLab 模型 可 以 模拟 植物 生长 过 程 中 的 生 
物 量 生 产 OG 和 需求 DD 的 变化 情况 。 四 个 玉米 品种 
灌溉 处 理 的 O 和 DD 值 均 大 于 雨 养 人 处理， 而 发 育 开 
始 时 灌溉 处 理 的 CD 比 均 大 于 雨 养 处 理 ， 但 叶片 
停止 发 育 后 ， 灌 溉 处 理 的 OD 比 均 小 于 雨 养 处 
理 (图 9)。 
3.4.4 玉米 植株 的 三 维 结构 模拟 结果 

基于 GreenLab 模型 可 以 根据 试验 测量 的 株 
高 数据 模拟 玉米 的 三 维 结构 。 图 10 给 出 了 四 个 
玉米 品种 两 个 处 理 的 收获 阶段 的 植物 三 维 结构 。 
由 图 可 见 ， 灌 溉 处 理 的 植株 高 度 均 高 于 不 灌溉 处 
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图 7 不 同 发 育 阶段 四 个 玉米 品种 两 个 处 理 的 叶 . 节 和 果 的 。 UA i 
GreenLab 5 DSSAT 模拟 结果 比较 EEE, HJ FY GreenLab 模型 计算 投影 面积 S, 


Fig. 7 Comparison of simulated leaves, internodes and fruits of 的 值 ， 并 与 试验 结果 一 致 。 Feng 等 C5: Fill FA 
different stage between DSSAT and GreenLab model for the four GreenLab 模 型 ， 采 用 差分 统计 的 方法 研究 了 玉米 
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e 间 具 有 良好 的 兼容 性 。Lemaire 等 ”基于 各 类 器 


ERLE utes mae nyo eats 单 株 到 群体 的 产量 转化 ， 计 算 结 果 与 群体 测量 结 
SEU E ow 果 一 致 。 这 些 研究 表明 ，GreenLab 模 型 仅 使 用 各 
s} ]ue "| 类 器 官 的 数据 ， 就 可 以 模拟 各 个 器 官 的 生长 发 育 
P jag og 动态 。 因 此 ，GreenLab 模 型 作为 基于 器 官 尺度 的 
e Jaa a 功能 结构 模型 ， 与 基于 过 程 的 作物 模型 完全 兼 
il j2 2 容 ， 弥 合 了 农艺 学 作物 模型 和 植物 学 构 型 模型 之 
DAN M E E. SOB 

叶 元 位 置 生长 周期 /GC 


Mi 2 本 研究 基于 DSSAT 的 模拟 数据 ， 利 用 

(a) 不 同 叶 元 位 置 的 叶子 、 节 和  (b) 不 同 生 长 阶段 叶子 、 节 、 果 和 m » "m" . 
果 的 重量 地 上 部 的 重量 GreenLab 模 型 可 以 很 好 地 模拟 出 LAI 和 地 上 部 总 
图 8 模拟 的 玉米 各 类 器 官 重 量 随 位 置 和 生长 阶段 的 变化 ” 重量 ， 与 试验 结果 一 致 。 此 外 ，GreenLab 模 型 计 
Fig. 8 Simulated variation of maize organ weights with phyto- 算 的 不 同 阶段 各 类 器 官 ( 叶片 : 节 、 果 ) 的 Äg r 
mere rand and different stages 量 也 与 DSSAT 的 模拟 数据 一 致 由 此 可 见 , DS- 
理 的 植株 ， 结 果 与 试验 测量 数据 一 致 GE), SAT 与 GreenLab 之 间 具 有 良好 的 兼容 性 。 初 步 
探讨 了 连接 传统 的 作物 模型 和 功能 结构 模型 的 方 
式 2 的 。 尽 管 不 一 定 改 进 基于 过 程 的 模型 的 模拟 
大 "m 结果 ， 基 于 器 官 尺度 的 模型 可 以 提供 关于 植物 结 
4.1 不 同 尺度 模型 连接 的 可 行 性 eee No 
构 的 更 深入 的 信息 77. DSSAT 中 的 测量 数据 是 


GreenLab 功能 结构 模型 与 基于 过 程 的 模型 之 各 类 器 官 的 信息 ， 而 在 GreenLab 中 ， 由 于 可 以 
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图 10 GreenLab 模拟 的 四 个 品种 两 个 处 理 的 玉米 三 维 结构 


Fig. 10 Simulated 3D structures of four maize cultivars with two treatments using the GreenLab model 


识别 絮 官 序列 ， 还 提供 了 对 单个 器 官 生 长 过 程 的 
模拟 (图 8 (a) ) 以 及 器 官 序列 的 逐个 周期 的 动 
态 变化 (图 8 (b) )。 此 外 ， 基 于 构建 的 Green- 
Lab 模 型 ， 可 以 进一步 进行 计算 实验 ， 优 化 源 库 
参数 ， 在 个 体 尺 度 上 研究 结构 对 产量 的 影响 "。 


4.2 不 同 尺度 模型 结合 的 应 用 领域 


虚拟 农业 中 的 一 个 研究 主题 是 直观 地 表达 施 
肥 灌 溉 等 措施 对 作物 个 体 生 长 的 影响 。 基 于 本 研 
究 的 从 群体 尺度 到 个 体 尺 度 的 研究 ， 可 以 基于 模 
型 快速 模拟 不 同 管理 措施 下 的 作物 形态 变化 ， 而 
无 需 逐 个 调整 器 官 大 小 。 

个 体 尺 度 的 植物 模型 往往 包含 器 官 尺度 的 子 
模型 ， 后 者 可 以 进一步 与 微观 尺度 的 遗传 模 
TY C9 或 者 光合 作用 模型 相 结合 2? 方便 地 引入 
单位 面积 的 光合 速率 对 个 体 生 长 的 影响 ， 从 而 充 
分 利用 作物 表 型 获取 的 信息 。 也 可 以 与 群体 尺度 
的 过 程 模型 (如 DSSAT 或 者 APSIM) 相 结 合 ， 
如 本 研究 所 示 。 个体 尺度 的 模型 可 以 作为 不 同 尺 
度 人 研究 的 桥梁 ， 进 行 更 广泛 程度 的 模型 研 
究 2 2。 应 用 领域 包括 气候 变化 情况 下 的 作物 生 
长 模拟 与 产量 产 期 预测 研究 ， 以 及 育种 研究 。 


5 结论 


本 研究 通过 基于 DSSAT 中 的 测量 数据 和 模 
拟 的 每 日 叶 、 节 和 果实 的 重量 数据 ， 利 用 Green- 


Lab 模 型 模拟 了 四 个 品种 两 个 处 理 的 玉米 植株 各 
个 器 官 的 生长 发 育 动 态 ， 探索 了 群体 尺度 DS- 
SAT 与 个 体 尺度 GreenLab 作物 模型 连接 的 可 能 
性 。 结 果 表 明 ， 基 于 共同 的 数据 结果 ,个 体 尺度 
的 GreenLab 模型 可 和 群体 尺度 的 模型 输出 一 致 
的 结果 。 除 了 可 复 现 模拟 叶 面 积 指数 和 地 上 部 总 
重 ， 还 可 进一步 反 演出 各 种 器 官 之 间 的 分 配 ， 从 
而 对 作物 的 生长 发 育 进行 更 详细 的 描述 ， 并 对 植 
物 的 生长 发 育 过 程 进 行 三 维 可 视 化 展示 ， 从 而 实 
现 从 群体 尺度 模型 到 个 体 尺 度 模型 的 转换 。 
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From Stand to Organ Level—A Trial of Connecting 
DSSAT and GreenLab Crop Model through Data 


WANG Xiujuan", KANG Mengzhen'*’, HUA Jing", Philippe DE REFFYE? 


(1. The State Key Laboratory of Management and Control for Complex Systems, Institute of Automation, Chinese Acade- 
my of Sciences, Beijing 100190, China; 2. School of Artificial Intelligence, University of Chinese Academy of Sciences, 
Beijing 100049, China; 3. Beijing Engineering Research Center of Intelligent Systems and Technology, Beijing 100190, 

China; 4. Qingdao Agri Tech Co. , Ltd. , Qingdao 266000, China; 5. CIRAD, UMR AMAP, Montpellier F-34398, France) 


Abstract: Crop models involve complex plant processes, which can be built in different scales of space and time, from molecule, 
cell, organ, tissue, individual to stand in space and from second to year in time. Based on different research requirements, 
switching the model scales can make the applicability of the model more extensive and flexible. How to switch the crop model 
from stand level to organ level is the content of this research. The DSSAT software (stand level) and functional-structural plant 
model 'GreenLab' (organ level) were chosen to explore the possibility to switch the crop model from stand to organ level. The 
DSSAT can simulate the growth and development processes of crops in detail according to the growth period by taking the data 
of weather, soil, crop management, and observational data as input. The GreenLab can simulate the growth and development 
and their interaction of crops by considering plant structure, and the model parameters can be estimated according to the mea- 
surements. In this study, the experimental data contains two parts: the measurements of four maize cultivars with two treatments 
(irrigated and rainfed) in DSSAT, and the simulations including the weights of leaves, internodes and fruits per day using DS- 
SAT based on the measurements. The simulation results of DSSAT were used to calibrate the parameters of the environmental 
(E), sink strength (P,), and remobilization (k, and k,) in GreenLab, and to compute the weights of leaves, internodes and fruits 
for each phytomer. The simulation results of the GreenLab model were compared and analyzed with the experimental data and 
the simulations of DSSAT. The consistency of calculation results could be an indicator to explore the method of building an in- 
terface between different-scale crop models, and to compare the characteristics of different models. The results showed that the 
GreenLab model could reproduce the simulation data of the DSSAT and the measurement data, including the leaf area index 
(LAD and the total weight of the plants, and further could compute the biomass for each organ (leaf, internode and fruit), and 
the biomass distribution among organs, the biomass production (Q), the demand (D) and the ratio between Q and D during the 
growth. Therefore, the detailed information of organ growth and development could be reproduced and the 3D structures of 
plant could be given. Finally, the advantages and application fields of different-scale model integration were discussed. 

Key words: crop model; different model scale; functional-structural plant model; model integration; DSSAT; GreenLab; parame- 


ter estimation 
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